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Tapered  Wiggler  Analysis  of  High  Grain 
Free  Electron  Laser  Oscillators 

I.  Introduction 

The  development  of  lasers  in  which  the  active  medium  is  a  relativistic  stream  of  free 
electrons  has  recently  evoked  much  interest.  The  prospect  of  operating  these  continuously 
tunable,  high  gain  devices  at  high  efficiencies  by  appropriately  tapering  the  wiggler  field 
parameters  has  significantly  increased  the  interest  in  these  devices.  The  dramatic  effect 
of  efficiency  enhancement  schemes  in  free  electron  lasers  relies  on  the  maintenance  of  the 
co-resonance  between  the  phase  velocity  of  the  ponderomotive  wave  with  the  decreasing 
velocity  of  the  electron  beam  as  it  gives  its  energy  to  the  radiation  field.  Variable  parameter 
wiggler  free  electron  lasers  have  been  previously  considered  for  operation  in  the  amplifier 
configuration  (Refs.  1-3)  and  also  in  the  micro-pulse,  low  gain  operation  in  the  oscillator 
configuration  (Refs.  4,5).  The  present  analysis  extends  the  range  of  application  to  the  case 
of  intense,  long  pulsed  electron  beams  in  which  the  self-consistently  derived  space  charge 
potential  (arising  from  the  axial  bunching  of  the  particles  in  the  beat  wave)  is  not  negligible. 
Also  included  in  this  analysis  are  arbitrarily  relativistic  beams  (no  ultra-relativistic  beam 
approximation).  Thusly,  the  simplifying  assumptions  that  7  >  1  and  (3.  1  are  not. 

incorporated  into  the  theory.  The  full  relativistic  effects  of  mildly  relativistic  beams,  which 
are  characteristic  of  beams  from  intense,  long  pulsed  devices,  are  adequately  considered. 
The  analysis  also  includes  the  high  gain  operation  in  the  oscillator  configuration  which  is 
characterized  by  exponential  field  growth  in  the  linear  regime. 

The  organization  of  the  paper  is  as  follows.  In  Sec.  II,  we  present  the  model  and  as¬ 
sumptions  used  in  the  numerical  simulation.  In  Sec.  Ill,  we  derive  approximate  t  heoretical 
expressions  for  the  nonlinear  efficiency.  In  .Sec.  IV,  we  present  a  summary  of  conclusions 
with  specific  examples  from  the  simulation  and  comparison  with  theoretical  expressions. 

II.  Simulation  Model 

We  have  conducted  a  fully  relativistic,  nonlinear  simulation  of  the  field  and  particle 
evolution  in  a  variable  parameter  wiggler,  high  gain  free  electron  laser  oscillator.  The 
simulation  model  which  is  kinetic  in  nature  is  coupled  with  Maxwell's  equations  to  self- 

consistently  include  the  evolution  of  the  space  charge  waves.  The  wiggler  field  we  consider 
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is  helically  polarized  in  the  plane  transverse  to  the  beam  propagation  direction,  and  is 
ideal  in  that  the  transverse  variations  in  the  vector  potential  are  neglected. 

Aw(z)  =  — au,(;)^  exp[-i  J  kw{z')dz')  c+  +  c.c.^,  (1) 

where  aw(z),kw(z)  are  the  arbitrarily  prescribed  wiggler  amplitude  and  wave  number  re¬ 
spectively  and  also  e±  =  (ex  ±  iiy)/2.  and  c.c.  denotes  complex  conjugate.  The  radiation 
field  model  is  also  one-dimensional  and  is  characterized  by  a  temporally  stationary  ampli¬ 
tude  with  arbitrary  spatial  dependence. 

As(z,t.)  -  <*,(2)exp[*(fc4z  -  e+  +  c.c.  (2) 

The  electrostatic  space  charge  potential  is  also  assumed  to  be  temporally  stationary  with 
arbitrary  spatial  dependence. 

${z,t)  --  01(c)exp[t{  /  ffc.  -I-  *„.( s' )|dc'  u>,f})  +  c.c.  (3) 

J  o 

The  vector  and  scalar  potentials  of  the  electromagnetic  wave  evolve  according  to  Maxwell’s 
equations  and  the  particles  evolve  in  accordance  with  the  kinetic  distribution  /  (Ref.  6). 
Transverse  effects  are  heuristically  included  by  considering  /  dr  J  dy  f 
oc  Trrg,  J  dx  J  dy  a,  ex  7ru>jj  with  rj,  equal  to  the  electron  beam  radius  and  w0  = 
(2 zn/kf)1/"  is  the  spot  size  for  an  oscillator  mode  of  Rayleigh  range  zp.  Inserting  the 
driving  currents  and  charge  densities  into  Maxwell’s  equations  and  assuming  i)  a  cold  elec¬ 
tron  beam,  ii)  vector  and  scalar  potentials  that  vary  slowly  as  compared  to  an  optical 
period  and  iii)  small  radiation  fields,  at  < <  au,  ,  one  obtains  the  following, 

dat  i  (  u>(,  \ 2  Ff3i0  .  ,  J  e~'v’  \ 

Tb  =2(7] 

d  i{k,  +  kw)\-  i/^bY  ih  0  /f~,0\ 

dz+  2  r  =  2  \  c  J  (k.  +  kw)\  &  /' 

where  we  have  introduced  the  dimensionless  potentials,  d.,,  a,,,,  <^,  corresponding  to 
|e|(a,,ou,,^i)/mc2.  The  nonrelativistic  beam  plasma  frequency  is  given  by  = 


Kit 
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(4ir»e2/m)1^2  and  F  =  t\/wq  is  the  filling  factor.  The  integration  over  initial  positions 
has  been  converted  to  an  integral  over  initial  phases,  i)>0,  with  {•  •  •)  =  (1/2tt)  J*w  ■  ■  •  <hj'0, 
and  the  phase  is  given  by  =  /„*  [fc,  +  kw(z')  -  vt/vz(z')\dz'  +  V’o- 

The  evolution  of  the  particle  phase  as  determined  by  the  Lorentz  force  equation  is 
given  by. 


»  dvz  _  ,  kt  dff~,  ( 1  ~b  fit  )k&kw  o w  , . 


il,"  _  jfc'  4.  r_i  _  i.»  _  _ _f  _ 

V  * +  t<2  dz  w  2/^h2  d- 


—  («,c  -  ase 


+  j*{£+W>  +  k»))+i  +  e-i*{£-i(k„  +  kw))ii  , 


where  prime  denotes  differentiation  with  respect  to  z.  The  oscillator  system  of  equations 
(Eqs.  (4,5))  are  solved  with  the  following  methods  and  boundary  conditions.  The  ordinary 
differential  equations  for  the  fields  and  particles  are  solved  with  a  four  point  Adams- 
Bashforth  predictor  corrector  scheme  which  is  initialized  with  a  fourth  order  Runge-Kutta 
algorithm.  After  propagating  an  arbitrarily  small  radiation  field  from  the  entrance  to  the 
exit  of  the  resonator,  a  portion  of  the  wave  (depending  on  the  reflectivity  of  the  cavity 
mirrors)  is  transmitted  out.  of  the  system  and  the  remaining  wave  is  used  to  initialize  the 
subsequent  pass.  In  this  manner  we  have  been  able  to  simulate  the  small  signal  exponential 
gain  to  the  nonlinear  saturated  evolution  of  a  tapered  wiggler  oscillator.  Representative 
results  are  presented  in  Sec.  IV. 


III.  Theoretical  Analysis 


The  simulations  are  aided  by  monitoring  a  diagnostic  quantity  which  may  be  obtained 
from  Poynting’s  theorem.  (Ref.  6) 


,  ^  _  Pr(z)  ~  Pr(0)  ,  , 

VR(Z)  =  - - =  Ve(  =  ) 


b(0))  -  b(=)) 

7o  -  1 


which  expresses  the  fact,  that  the  energy  loss  by  the  particles  must  appear  in  the  radiation 
field.  These  quantities  are  monitored  at.  every  axial  grid  of  the  simulation  and  provides  an 
excellent  check  on  the  consistency  of  the  field  and  particle  evolution. 
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Defining  the  saturated  state  as  the  condition  that  the  radiation  power  gained  per  pass 
equals  the  radiation  power  loss  at  the  mirrors,  one  obtains 


I  0.(1)  I1  -  1 0.(0)  |!  ]  =  (1  - 


.(£) 


or 


PR»ai(L) 


VR(L)Ph 

(1  -R)' 


(8) 


where  PRtat  is  the  radiation  power  at  saturation,  Pi,  is  the  input  electron  beam  power,  tjr 

■i 

is  the  single  pass  radiation  efficiency  as  defined  in  Poynting’s  theorem  and  R,  is  the  power 
reflection  coefficient  for  the  cavity  mirrors.  This  expression  aptly  indicates  the  possibility 
of  achieving  high  radiation  powers  in  high  Q  resonators  with  modest  beam  powers.  The 
remaining  theoretical  sections  are  devoted  to  estimating,  t)r(L),  the  single  pass  radiation 
efficiency  in  terms  of  the  prescribed  beam  and  wiggler  parameters. 

The  evolution  of  the  radiation  power  can  be  obtained  from  the  vector  potential  equa¬ 
tion  (4), 


~dT+t<t>-  |a'  l=i(7} 


(9) 


where  <t>.  is  the  radiation  phase,  a$  =|  a,  \  c**’  and  0,  —  V’  +  4>>-  Multiplying  the  real  part 
of  Eq.  (9)  by  |  a,  |  (c  j  Air)k2iriv1  one  obtains 


dPR(z) 

dz 


A 


7o  -  1 


k$d-i 


(10) 


III.  a)  Resonant  Particle  Approximation 


The  simplest  estimate  of  the  single  pass  efficiency  as  considered  by  Refs.  1 .2  can  be 
obtained  by  making  the  Compton  and  resonant  phase  approximations.  The  resonant  phase 
approximation  of  Eq.  (10)  is  given  by, 


dPn(z) 

dz 


A 


7o  -  1 


k$d i 


sin  iI',r 


(ID 
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y***M'CV *  "r 


*'t  l 


The  sine  of  the  resonant  phase  is  defined  by  evaluating  the  fixed  point  of  the  pendulum 
equation  (Eq.  (6))  at  the  resonant  energy  given  by  7^  =  k,fi2 /f3z(l  4-  /3z)kw,  with  /r  = 
1  4-  a2w. 


0  =  <t>"  +  k'w  - 


da2w 


(1  +  fiz)kwks . 

a, 


a,  |  sinV’aR. 


(12) 


2 PWn  dz 

Neglecting  the  rate  of  change  of  the  radiation  phase  as  compared  to  the  wiggler  period 
and  inserting  the  results  from  Eq.  (12)  into  Eq.  (11)  yields, 


dPR(z) 

dz 


=  Prto(:h 


(13) 


where  tj'0  =  [(if/rtfr - 2^  Equation  (13)  can  be  integrated  analytically  and 

yields  the  following  efficiency  for  the  resonant,  particle  approximation, 


V Re,  =  /  Vo (z)dz  - 
Jo 


Pr(L)  -  PR( 0)  _  2 13l0  7r(0)  -  7 r(L) 


(14) 


Pb  (1  +  fizo)  "70-1 

This  efficiency  is  the  relativistic  generalization  of  the  bucket  efficiency  of  Refs.  1,2  and 
yields  efficiencies  which  are  of  the  order  of  70  percent  larger  than  the  efficiencies  obtained 
from  the  simulation  of  continuous  electron  beams.  Although,  when  the  simulations  are 
performed  for  prebuncliecl  beams  characterized  by  macro-particles  injected  at.  the  resonant 
phase,  Eq.  (14)  yields  very  good  agreement,  with  that  simulation  efficiency. 


III.  b)  Deeply  Trapped  Particle  Compton  Limit 

In  order  to  better  approximate  the  ensemble  average  over  initial  phases,  the  assump¬ 
tion  of  a  macro-particle  at  the  resonant  phase  will  be  replaced  with  an  explicit  average 
over  the  orbits  of  deeply  trapped  particles.  As  a  first  example  we  shall  consider  again  the 
Compton  approximation.  Expanding  the  pendulum  equation,  Eq.  (6),  about  the  resonant 
phase,  =  4'tR  +  H'  at  the  resonant  energy  7^  =  fi2kji3z(l  -!  liz)kw  yields, 


W  =  <4  -  -  V,»r( ~ )  -  ~ r'"— I  I  {  sin  V’,R  +  ros  V’,R  *</’},  ( IS) 


2 PIUr  d: 


KTr 
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where  we  shall  explicitly  retain  the  drift  and  acceleration  of  the  resonant  phase,  ^'aR  ^ 
0,  4,"r  #  0.  Using  the  definition  of  the  sine  of  the  resonant  phase,  Eq.  (12),  the  pendulum 
equation  becomes, 


H"  =  -4\r{~)  -  k;yn{z)6il\  (16) 

where  k2ayn(z)  =  (1  +  /3z)kwkaaw  |  aa  |  cos  4'>r/PIir-  The  oscillation  of  the  deeply  trapped 
particle  can  be  explicitly  solved  by  multiple  scale  techniques  which  exploit  the  desperate 
time  scales  of  the  synchrotron  bouncing  of  the  particle  and  the  adiabatic  variation  of 
wiggler  parameters.  The  solution  is  given  by, 

*4’ = +  £(4>o,z)cos  f  kayn(z')dz' -  p(z)  sin  /  kayn(z')dz'  (17) 
'Synt-J  Jo  Jo 

where  £(V’o,*)  =  y/k,y„(0){4'o  -  4'*r( 0)  +  4'aR(0)/k2ayn(0))  /  y/kay„(z),  p(z)  = 

V’ifl(O)/ \A’»y*»(0)^*yn(~)  and  we  have  chosen  boundary  conditions  corresponding  to 
V’»(0)  =  V’o  and  V’1(0)  =  0.  Making  use  of  this  solution  for  the  ensemble  average  yields, 


k~  —  oo 

e-« (18) 


f=- 


where  Jm(x)  is  the  Bessel  function  of  the  first  kind  of  order  m.  By  neglecting  the  oscilla¬ 
tions  at  the  sideband  harmonics  one  obtains, 


(sin^4)  =  sin (tf.n(s)  -  4'"R(:)/k2yn{z))J0(p(z))^,  (19) 

with  $  =  j  -~Jo{S(4'o,z)), 

J  rfV70(V)  =  A70(A)  +  y  [F0(A)7,(A)- F1(A)J0(A)], 


and 


where  Hk{%)  is  the  Struve  function  of  order  k.  The  results  of  the  ensemble  average,  when 
included  in  the  power  evolution  equation,  yields, 


- "T^  ~  =  VlPb  -  Vi'nVhPRiz), 

az 

where  rf,  =  3  J0  (p(z))  cos  % ( ~ ) 

^  a  y  n  *  ~  * 


(20) 


and 


7  R 


Vi i  =  —rr^Jo  (*>(  =  ))  sin  pr TT,  cos 

70  t  ynl-/  1  C.R 


,  ,  .««■“>&  (  7u  "  1  \ 

Vsr(~) — —  -TTr 

n«  V  fltF  / 


1/2 


Although  Eq.  (20)  can  be  solved  analytically  for  the  case  of  constant  i]'r  and  the 
transcendental  solution  for  Pr(z)  makes  it  impossible  to  evaluate  the  single  pass  efficiency, 
i]r(L)  -  (J°!r(0)  PB(L))/Ph.  For  this  reason  we  choose  to  retain  the  slow  variation  of 
r)'j  and  i]'j i  and  solve  Eq.  (20)  numerically  and  compare  the  results  to  the  simulations. 
This  efficiency  yields  results  that,  are  within  10  percent  of  the  simulation  efficiency  for  the 
Compton  regime  operation.  For  parameters  of  physical  interest.  ij'f  "•>  rf'n  and  the  evolution 
of  Eq.  (20)  is  dominated  by  the  first  term.  The  deeply  trapped  particle  approximation  to 
the  efficiency  is  dominated  by  the  bucket  efficiency  times  the  product  of  three  form  factors, 
3 ,  Jo( p(  ~ )  land  cos  V’"R(  r  )/k2ty„ ( s ),  which  are  due  to  the  average  motion  in  the  bucket  ,  drift 
and  acceleration  of  the  bucket  respectively. 

III.  c)  Deeply  Trapped  Particle  Raman  Limit 

The  space  charge  effects  can  be  easily  included  with  this  technique  of  deeply  trapped 
orbit  averaging.  The  pendulum  equation  in  the  presence  of  space  charge  waves  and  ex¬ 
panded  about,  the  resonant  phase  and  energy  is  given  by. 


where  p  is  the  phase  of  the  electrostatic  potential,  and  we  have  neglected 

the  rate  of  change  of  the  space  charge  potential  as  compared  to  the  ponderomotive  wave¬ 
length.  By  making  use  of  the  resonant  phase  definition  in  the  Compton  regime,  k'u,  - 
^3 /( 2/3? 7 p  )da*„ / dz  =  (1  +  /?;)fcu,fcsdu,|a4|  sin  4'sr/ the  pendulum  equation  can  be 
cast  in  the  form, 

S4"  =  a(z)  ~  k2  J (22) 


where  a(:)  =  -4"R  4- sin(V’,R  4  <p  <p„), 

"zTr 


7.2  2ks{\  /T  )  .  .  i  I  ,  ,  /  4  |  (  4-  rfz  )kwk.,  . 

"'syn\  "  )  -  "13  2 - 1*  ’  +  «-'«•) (<Pl  I  c°s(  4\,R  f  -  <f>,  )  H - rj— - «  |o,  |  COS  i/\,R  - 

Pz~)~R 

The  magnitude  of  the  electrostatic  potential.  |<£]|,  and  the  phase  shift  between  the 
electrostatic  and  radiation  phases,  y?  -  </>,  can  be  estimated  from  the  evolution  equation 
for  the  space  charge  potential. 


I  4>\  1=  —  (cos(V\,  +  y?  -  4>.<)  i  sin (</’.«  4  s?  </>,))•  (23) 

(«3  +  ^«.)‘<"' 

Since  the  left-hand  side  of  Eq.  (23)  is  real,  one  has  that  (sin(  t/\«  +  s?  -  </>.« ))  --  0.  Using 
the  deeply  trapped  orbits  of  the  Compton  approximation  to  evaluate  the  ensemble  average 
yields, 

(sin(V’s  +<p  -  <£,))  ~  (p(  ?))  sin( «/’.,/}  4  y?  <t>,  v"B/k;  )  0.  (24) 


V’.iR  +  ?  -  *• 

Thusly,  the  coefficients  in  the  space  charge  pendulum  ecjuation  are  given  l>y. 


(26) 


k‘‘niz)  ~  V  +  k,)c-  *J° (,,(l)l  “s 


,  (1  +0: )ktkw  .  , 

+  - mT~2 - att)|a»l  COS 

PzTr 


where  the  space  charge  potential  modifies  the  synchrotron  wave  vector  k,yn  by  increas¬ 
ing  the  synchrotron  period  due  to  electrostatic  repulsion,  and  also  increases  the  effective 
bucket  acceleration  a(z).  Again  making  use  of  multiple  scale  techniques  the  solution  to 
the  pendulum  equation  is  given  by, 


**  =  j?r~+i(4'9,z)coS  f  k„yn( ~  )d,z'  —  p{z) sin  f  ksyn(z')d.z ',  (27) 

JO  JO 


where  £(V’o,=)  =  yjk, yn(0)(^0  -  4'»r{ 0)  -  a(0)/fc;yn(0))/y  £,,,„.(.:)  and  p(r)  = 

#*(  0)/y/i  ayn  (0)fcsyn(2).  Performing  the  ensemble  average  with  these  orbits,  the  power 
evolution  equation  becomes, 


=  spViPb  +  sPv'nVP>>PR(-h 


where  =  3J0  (p(r))cos 

3 yn\~ ) 


,  #  7 fi  u  .  «(-)  .  awu>t  ( 70  "  1 Y 

and  Jo  (/>(-))  sin  cos  0jR  -T-  (  -r~—  ) 

7-1  ^n(.)  aR  V  ) 


and  3  =  J  —■ Jo  (£(4'o ,  c ))  . 


The  numerical  solution  to  Eq.  (28)  is  compared  to  the  simulation  results  in  the  space  charge 
dominated  regime.  The  comparison  presented  in  Sec.  IV  is  characterized  by  agreement 
between  theory  and  simulation  on  the  order  of  10  percent. 


IV.  Results  and  Summary 


The  simulation  model  is  capable  of  tracking  the  evolution  of  the  radiation  field  from 
the  small  signal  exponential  gain  through  the  nonlinear  saturated  state  of  the  oscillator. 
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During  the  small  signal  evolution  of  the  fields  for  which  the  growth  rate  is  constant,  the 
growth  rate  from  the  simulation  is  compared  to  the  growth  rate  obtained  from  the  solution 
of  the  linear  dispersion  relation.  This  comparison  is  presented  in  Fig.  1  for  the  case  of, 
beam  current  I  =  300.4,  beam  energy  7  =  2.5,  wiggler  field  strength  Bw  =  1000G',  wiggler 
period  (w  =  4.0cm,  filling  factor  F  =  0.2  and  no  taper  of  the  wiggler  amplitude  or  period. 
The  comparison  of  the  simulation  and  theory  yields  excellent  agreement  throughout  the 
unstable  spectrum. 

An  example  of  the  small  signal  to  nonlinear  evolution  of  the  fields  is  presented  in  Fig. 
2.  In  this  example  the  field  evolves  over  200  radiation  bounce  times  for  the  case  of  beam 
current  I  =  3004,  beam  energy  7  —  2.5,  wiggler  field  strength  Bw  =  1000G',  initial  wiggler 
period  (w  =  4.0cm,  filling  factor  F  =  0.2,  linear  taper  of  the  wiggler  period  of  50  percent 
and  mirror  reflectivity  R  =  0.98.  The  saturated  power  and  efficiency  are  1 ,9G'IF  and  17 
percent  respectively;  and  exhibits  the  scaling  indicated  in  Eq.  (28). 

The  single  pass  efficiency  exhibits  a  weak  dependence  on  the  radiation  power  at  the 
entrance  of  the  resonator.  In  Fig.  3  we  present  a  comparison  of  the  single  pass  efficiency 
from  the  simulations  and  theory  as  a  function  of  the  input  power.  As  noted  in  Sec. 
III.  b),  the  theoretical  evaluation  of  the  efficiency  neglects  oscillations  at  the  harmonics 
of  the  synchrotron  period.  The  simulation  model  contains  the  complete  dynamics  of  the 
electrons  as  dictated  by  the  pendulum  equation  and  the  simulation  efficiency  tliusly  exhibit  s 
oscillations  at  synchrotron  harmonics  with  an  amplitude  of  the  order  of  the  well  depth. 
The  bars  on  the  data  points  from  the  simulation  indicate  the  amplitude  of  this  oscillation 
in  the  efficiency.  The  resulting  comparison  indicates  agreement  on  the  order  of  10  percent 
between  the  simulation  and  theoretical  expressions. 
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